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Abstract 

Single top quark production in electron photon interactions provides 
a clean environment for the measurement of the Cabbibo-Kobayashi- 
Maskawa matrix element Vtb- Aiming an experimental precision at the 
percent level the knowledge of radiative corrections is important. In this 
paper we present results for the radiative corrections in quantum chromo- 
dynamics. 



1. Introduction 



One of the fundamental unsolved problems of todays high-energy physics is the exact mech- 
anism of electroweak symmetry breaking (EWSB). Due to the fact that only the top quark 
couples with a Yukawa coupling of order one to the so far unobserved Higgs boson it is natu- 
ral to assume that the top quark plays a special role in the context of the EWSB. For example 
in so-called dynamical symmetry breaking models the scalar Higgs field — responsible for 
the spontaneous symmetry breaking in the Standard Model — is replaced by a composite 
scalar operator. Such an operator could be built for example from heavy fermion fields. 
Having eliminated the elementary scalar field from the theory the problem of the large mass 
corrections due to quantum corrections is solved. Examples for such models are techni- 
colour models (for a review see ref. Q and references therein), top-condensate models EL 
and top-colour models BHU. For the search of such extensions a precise understanding of 
the top quark sector of the standard model is necessary. 

At hadron colliders both single top quark production as well as top quark pair production 
have been studied extensively in the past. The differential cross section for top quark pair 
production is known to next-to-leading order (NLO) accuracy in quantum chromodynamics 
(QCD) 13 [3 03 . In addition the resummation of logarithmic enhanced contributions has 
been studied in detail in refs. JHEIEIEIE 1- Recently also the spin correlations between 
top quark and anti-top quark were calculated at NLO in QCD ifLUl . Due to the fact that single 
top quark production provides an excellent opportunity to test the charged-current weak- 
interaction of the top quark it has also attracted a lot of interest in the past. In particular 
NLO corrections were studied in refs. I1T51 [Tol [P71 ITSll. In ref. [18| the NLO corrections 
for the fully differential cross section is given keeping also the spin information of the top 
quark. On the experimental side the situation is not very conclusive for the moment as far 
as single top quark production is concerned. Due to limited statistics at run I of the Tevatron 
collider only upper bounds were obained in ref. [19]. In ref. ll20ll the possibility to measure 
the electroweak couplings in single top quark production at the LHC is studied. In particular 
also the sensitivity to new physics is discussed. 

As far as lepton colliders are concerned much effort has been devoted to top quark pair 
production in e + e -annihilation. In particular the total cross section in the threshold region 
is known at next-to-next-leading order (NNLO) in QCD. (For an overview of the theoretical 
status we refer to ref. ED-) Momentum ll22ll23ll as well as angular ll24l 1251 distributions 
were studied in detail. In the continuum region the total cross section for massive quarks is 
known to order oc,$ in the coupling constant of the strong interaction [26|. In order a,, 3 the 
quartic mass corrections to the total cross section are also known ll27ll . Also less inclusive 
observables have been studied in great detail. For example the spin structure of top anti-top 
system is completely known at order a s |EH1 1211110111111121 and partially known at order a s 2 
fl33l . Futhermore the 3-jet results obtained for massive ^-quarks Il34ll35ll36ll37ll38ll39l are 
also applicable to top quark physics 11401 . 



1 



Less attention has been devoted to single top quark production. Studies at tree level can be 
found for example in refs. EH EH EH EH Ell. In the refs. EHUD special emphasis was put 
on single top quark production at LEP II. For a top quark mass arround 175 GeV the pro- 
duction rates in the standard model are to small to be detectable at LEP II. This was recently 
confirmed by the L3 collaboration [46 1. In ref. ll45l also the single top quark production in 
electron photon collisions is studied. The electron photon reaction provides a clean environ- 
ment for the study of single top quark production because there is no background from top 
quark pair production. As a consequence this reaction is very well suited for the measure- 
ment of the weak couplings of the top quark. In particular it has been shown in ref. P31l that 
using polarised beams the Cabibbo-Kobayashi-Maskawa (CKM) matrix element V t b can be 
measured with an uncertainty of 1% at the 2a level. In this analysis 10 4 top quark events 
were assumed which corresponds to a luminosity of 100 fb 1 of an electron photon collider 
operating at y/s = 500 GeV. 

Aiming an accuracy at the percent level it is clear that the knowledge of the QCD corrections 
is mandatory. This is the main purpose of the present paper. In addition we study the struc- 
ture of logarithmic enhanced contributions which are related to initial state singularities. The 
full dependence on the b-quark mass is kept. This allows a systematic comparison between 
the structure function approach and the fixed-order calculation. Furthermore close to the 
threshold effects of the finite £>-quark mass are important. To calculate the QCD corrections 
we use the effective W-approximation. In the VK-approximation the scattering process which 
needs to be studied is W + y — > tb. Using the ^-approximation which describes the momen- 
tum distribution of the W-boson in the electron a prediction for e + y — > tb\ e can be obtained. 
A more detailed discussion is given in section |5] 

The outline of the paper is as follows. In section[2]we discuss the kinematics and the leading 
order results for the reaction W + y — > tb. The virtual corrections to this process are discussed 
in section |3 In section |4] the real corrections are calculated. In particular the cancellation 
of the infrared singularities is shown. In the following section we present the results for the 
subreaction W + y — > tb as well as for the reaction e + y — > tbv e . We finally close with our 
conclusions in section[6] 



2. Kinematics and leading-order results 
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Figure 2.1: Feynman diagrams for W + y — > tb in leading-order 
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In the following we study the reaction 

W + (k w ) +y(k y ) -> t(k t ) +b(k b ), (2.1) 

where we treat both outgoing quarks as massive. For later use it is convenient to define 
dimensionless variables. In particular we define the rescaled masses 

« = "t (2.2) 
s 

and the energy fractions 

2(k-ki) 

Xi = -± x l (2.3) 

s 

with k = k w + ky and s = k 2 . For the reaction given in eq. (12.11) the energy fractions are fixed 
completely by the kinematics: 

x t = l+Zt — Zb, and x b = l+zb-z t . (2.4) 

This is no longer true when the emission of an additional gluon is considered (c.f. section 
El). The calculation of the Born matrix elements for the reaction (12.11 ) is straightforward so 
we just quote the results here. The corresponding Feynman diagrams are shown in fig. l2.il 
For the W boson we distinguish between longitudinal and transverse polarization while for 
the photon we average over the incoming polarization. For partons in the final state the 
polarization (and the colour) is summed over. In terms of the leading-order squared matrix 
elements \M Q ' | for transversely /longitudinally polarized VK-bosons the differential cross 
sections are given by 

da TL _ NX bt 1 v \<u T > L \ 2 

where X bt is defined as 

X bt = X(l,Zb,z t ) (2.6) 

with 

X(x,y,z) = \/ x 1 +y 2 + z 2 - 2xy - 2xz - 2yz. (2.7) 

The number of colours is denoted by N. For transversely polarized W-bosons the normal- 
ization 5\£r is given by 9\£r = 2-2. For longitudinally polarized VK-bosons we have fVz, = 2. 
The squared matrix element for longitudinally polarized W-bosons is given by 

V- i„,L|2 2k f &(zb 3 -2zb 2 z t + ZbZ 2 + (3zb 2 -ZbZ t )z w ) 

L \ M o I ^ 

y,t,b Pol. 



9z w (l — z w ) 2 V (x b — X bt z) 2 



32(zb 2 z t -2z b zt 2 +zt 3 + (3zt 2 -z b zt)z w ) l\ l\ 



(x t + XbtZ) 2 x b~ ^btZ X t + XbtZ 

+ l2 z X ht (3zb - 3z t -z w )z w + lSz 2 X 2 bt z w 2 + 18z h 2 - 9z b [l + 4z t ) - 9z t 
+ lSz 2 + z w (36z b 2 + 12%(1 - 6z t ) + 60z t + 36z 2 ) 

+ z w 2 (2-l8zb 2 + 39z t -lSz t 2 + 3zb(5 + l2z t ))), (2.8) 
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where z denotes the cosine of the scattering angle 

Z = COS Qby = — cos f y (2.9) 

in the center of mass system and the prefactor K is given by 

K=sV2naGf\V tb \ 2 z w s. (2.10) 

Here G/ denotes the Fermi constant and V tb the Cabbibo-Kobayashi-Maskawa (CKM) ma- 
trix element. The functions l\, l\ are given by 

l\ = 2(-\0z b 3 -2z b 2 (l-Uz t )+z h (\+4z t -Uz t 2 )+Zt-2z t z + 2 Zt 3 

-2 Zw {\9z b 2 + z t 2 + A Zb ( 1 - 3z t ) ) + zj(zt -7z b )), (2.11) 



%{z b - 2z b 2 + 2 Zb 3 + (1 +4z b - Sz h 2 )zt - 2(1 - 5z b )z t 2 - A Zt 3 
+z w {-2 Zb 2 - 2(4 - 9z b )zt - 20 Zt 2 ) +z w 2 (z b - 7z,)). 



(2.12) 



The above matrix element is given in 4 dimension. In the context of the QCD corrections we 
will also need the squared matrix elements in d dimensions. The matrix element in d = 4 — 2e 
dimensions for longitudinally polarized W-bosons is given by 



E 

y.t.b Pol. 



94, 



L |2 
0,dl 



2(zb + zt) 



y.t.b Pol. 



i4 



+ 



8 



9 \x b -XbtZ x t + XbtZ 



In the derivation of the results above we used 

e1(k w )el(k w y- 



KV k w % v + k w v ky» m w 2 kfkj 



(k w ■ k y ) 



(k w -k y ) 7 



(2.13) 



(2.14) 



with £l(&w) being the polarization vector of the incoming longitudinally polarized W boson. 
The squared matrix element for transversely polarized W -bosons reads: 



E 

y,t,b Pol. 



+ 



2k / 64z f z w (zb-Zt-Zw) 1 6zbZ w (zb-z t +z w ) 



9z w (l-z w ) 2 \ {x t +Xbtzy 



(xb-lbtzY 



+ 



l T 



Xb — XbtZ x t + XbtZ 



■ + 6zl bt (z w (l - 3z b + 3z t ) + 6z w 2 (z t - Zb) 



+ z w 3 ( 1 + 3z h - 3z t ) ) - 9z 2 X 2 bt z w ( 1 + z w l ) + z w ( - 1 9 - 45z h z - 1 Sz t 



A5z t 2 + z b ( 18 + 90z f ) ) + z w 2 (36z b 2 -z b (l2 + 12 Zt ) - 60z t + 36z t 2 ) 



+ z w 3 (-l9- 9z h 2 + 6z t - 9z t 2 -Zb(6-lSzt)) 



(2.15) 



with 



/[ = 4z w (l + Uzb 2 + 2zb(l-&z t )-2zt + 2z t 2 + 2z w pZb-z t )+z w 2 ), (2.16) 



l T 2 = l6z w {l+2zb 2 -Zb(2 + 10z t )+2z t + 8z 2 -2z w (zb-4z t )+z w 2 ), (2.17) 
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k w % v + k w v k/ m w 2 k/kj 



where we have used 
for the sum over the two different transverse polarizations of the W -boson. 



(2.18) 



For later use we also give the squared matrix element for transversely polarized VK-bosons 
in d dimensions: 



I 

y,t,b Pol. 
1 



T\2 



= E 

y,t,b Pol. 



T\2 



+ £K 



2z% t --zht(l + Xzb-zt)) 



( 32zb(zb + zt-z w ) 128zt(zb + z t -z w ) 



9(l-z w ) 2 \ (x b -htz) 2 



(Xt+hbtZ) 2 



+ 



16 



(1 - Zb - 5z h 2 - ( 1 + 4z b )z t + zt 2 -(l-5z h +zt)z w + z w 2 ) 



Xb - ^btZ 

+ — -» — ( 1 - z h + z b 2 - ( 1 + Zb)z t - 2z 2 -(l+z h -2zt)z w + z w 2 ) 
x t + HtZ 

- 2 (-36(zb + z t )+31(l + z w 2 )-3Sz w + 3(l - z w ) 2 {2 + 3 ( Zb - z t ))(z b - z t )) 



- 4K£ 



1-1 



+ ■ 
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9 \x b -'kbtz Xt + XbtZ 



(2.19) 



Performing the remaining integration over the scattering angle we obtain the leading-order 

T L 

total cross section c t0 ' t for transversely/longitudinally polarized W-bosons: 



t,l_ 1 aG f \V tb \ 2 N 1 / r ,L 



'tot 



9\[2 (1-Zw) 3 ^r,L 



(2.20) 



with 



K L = 2l bt (28zb 2 -9z t + 28z 2 -Zb(9 + 56z t )+2z w (l%Zb 2 + 41z t + 18z 2 

+z b (5 - 36z t )) - Zw 2 (-8 + I2z h 2 - 3z fe (l + 8z f ) - 27z f + 12z, 2 )) (2.21) 
*T r = 4z w X fe ( - 2A Zb 2 + 6z b (3 + Sz t ) - 1 1 - I2z t - 2A Zt 2 - 2 Zw {-9z b 2 + 3z b (l + 6z t ) 
+5 + \5z t - 9z t 2 ) - z w 2 (6z b 2 - \2z b zt + 1 1 - 6z t + 6z t 2 )) . (2.22) 
The velocities (3 f , % of the outgoing quarks in the center-of-mass system are given by 

^bt , r> ^bt 



Pr 



and 13^ 



(2.23) 



l+Zt-Zb 1+Zb-Zt 

Note that the relation between the cross sections for unpolarized, transversely and longitudi- 
nally polarized W-bosons is given by 



<5= 1 -{2<5 T +<5 L ). 



(2.24) 



Furthermore we note that the structure of the logarithmic terms in eq. (12.201) is universal and 
can be obtained without an explicit calculation. In particular the singular contribution in the 
limit m b — ► can be written as 

m -»0 f l 

a(W + (k w )y(ky) -> t(k t )b(k h )) / dxf bh (x,m b 2 ,s) x o{W + (k w )b{xky) -> f&)), 

(2.25) 

where f b /y(x,m,b 2 ,s) can be interpreted as the bottom distribution in the photon (at scale s): 

f b/y (x,m b 2 ,s) = ^^Wln f^J (2.26) 

with the Altarelli-Parisi Py^ q q{x) kernel [47 1 given by 

P y ^ qq (x)=x 2 + (l-x) 2 . (2.27) 

(We denote with Q b the electric charge of the b-quark in units of the elementary charge e.) 
A more detailed discussion will be given in section |5] where the so-called structure function 
approach is investigated. 



3. Virtual corrections 

In this section we discuss the calculation of the virtual corrections. In particular we sketch 
briefly a few technicalities of the calculation, discuss the ultraviolet (UV) and infrared sin- 
gularities (IR), and carry out the renormalization. 

We work in renormalized perturbation theory, which means that the bare quantities (fields 
and couplings) are expressed in terms of renormalized quantities. By this procedure one ob- 
tains two contributions: one is the original Lagrangian but now in terms of the renormalized 
quantities, the second contribution are the so-called counter terms: 

L^o.Ao.mo.go) = L( x P R ,A R ,m R ,g R )+L ct ( } ¥ R ,A R ,m R ,g R ). (3.1) 

The first contribution yields the same Feynman rules as the bare Lagrangian but with the 
bare quantities replaced by renormalized ones. In the following we renormalize the quark 
field and the quark mass in the on-shell scheme. The conversion of the on-shell mass to the 
frequently used MS mass or to any other renormalization scheme can be performed at the 
end of the calculation. In spite of the fact that the calculation presented here is a one-loop 
calculation, it is still leading-order in the strong coupling constant a s . As a consequence the 
renormalization of the coupling constant does not appear. 

Whereas in the MS scheme the renormalization constants contain only UV singularities, in 
the on-shell scheme they contain also infrared divergences. We use dimensional regular- 
ization to treat both types of divergencies. Although at the very end all the divergences 
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must cancel it is worthwhile to distinguish between UV and IR singularities so that one can 
check the UV finiteness after renormalization and the cancellation of the IR singularities 
independently. To start with let us discuss the contribution of the one-loop diagrams before 
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Figure 3.1: Virtual corrections to W + j — > tb 



renormalization, that is the contribution from L (^r, . . . ,mR,gR). The corresponding Feyn- 
man diagrams are shown in fig. 13.11 The calculation of the one-loop amplitude is tedious but 
straightforward. To reduce the one-loop tensor integrals to scalar one-loop integrals we used 
the Passarino-Veltman reduction procedure ll4*8l . The two one-loop box integrals are given 
by 

D d = J_ f d d ( \ 

0)1 in 2 J (£ 2 -m t 2 )(£-k t ) 2 ((£-k t -k b ) 2 -m b 2 )((£-k t -k h + kj) 2 -m b 2 y 



u 0,2 — u 0,l 



(3.2) 

(m t ,k t )^(m h ,k b ) 



All the simpler topologies follow from these integrals by dropping one, two or three propa- 
gators. For example, using the standard notation [48 1 the infrared divergent triangle integral 
is given by 

If 1 

Co (1,2,3) = J d £ ^2 _ ^_ kt) 2^_ kt _ kb) 2_ mb2) ■ 0.3) 

The one-point, two-point and the finite three-point integrals have been calculated in the stan- 
dard way using Feynman parameterization. We have checked that our results for these inte- 
grals agree with the numerical evaluation given by the FF package of G. J. van Oldenborgh 
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Il4"9ll50l . Here we give only explicit results for the two IR divergent integrals. The triangle 
integral is given by 



(27i//) 2£ ReCg(l,2,3) = -^r(l+e)(^l -ln(p) - -21n(X fa )ln(p) 



1 / Ann 



shit ' V s 



,2\ e 



1 



£ vr/ 3 



^-(ln((o f ) 2 + ln((o fc ) 2 ) - Li 2 ((O b ) - Li 2 ((O f )] + 0(e) 

l . [A%n 2 \ z l 



r(l+e) — ^ -ln(p)+C (l,2,3), (3.4) 



with 



and 



sXtt V s I £ 



i — Zt - Zb + l^bt 

To calculate the infrared divergent box integral we have used two different methods. First we 
have considered a subtracted version of the integral which can be calculated in 4 dimensions 
using Feynman parameterization. From this result the desired result for the box integral can 
be easily obtained. The second method is based on the result given in ref. iBTl . There the 
infrared singularity is regulated by a small photon mass X. This result can be converted to 
dimensional regularization using the substitution 11521 

ln(X 2 ) -> i - y+\n(4n/j 2 ). (3.7) 
We found agreement of the results obtained by the two methods. The final expression reads 

s J 



x 



lln(p )+ ln(p)ln( ^-^ )-ln(x3^ 



+21n {j^) ln (*3) -2Li 2 (-px 3 ) +Li 2 (p 2 ) - 2Li 2 (-p/x 3 ) + ^ 2 }, (3. 



8) 



with 



s i . -w w Zw Zt Zb XyZwi Zt ^ Zb) /o n\ 

*3 = ~z {Z w -Z t -Zb + Hz w ,Zb,Zt)) = \ z~, r- (3.9) 

2m b m t y z w - z t - Zb - A(z w , z t , z b ) 

In the actual calculation we have replaced the box integral in d = 4 — 2£ dimensions by the 
box integral in (d + 2) dimensions. This can be done by the use of the relation 

Di 7 = -^D d + 2 (3.10) 
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where D^ 7 is the coefficient of the metric tensor g/j V in the decomposition of the four-point 
tensor integral D'L (cf. eq. (F.3) in ref. EU), which in turn can be expressed as a linear 
combination of the scalar box integral Dq and scalar triangle integrals Cq (i, j, k) in d dimen- 
sions. This procedure has the advantage that only one infrared divergent integral appears. 
As a consequence the extraction of the IR divergent contribution is simplified. Very often 
this procedure yields also a reduction of the algebraic complexity of the coefficients of the 
scalar integrals. Using Dq +1 instead of Dq we obtain the following result for the contribution 
involving Cq(1,2,3): 



si^l 2 



= ^C F (z, + z,-l)(2^) 2e ReC (l,2,3)^ £ \M Q T ^\ 2 (3.11) 

IR " div - n y,t,bPoL 



(In the following discussion of the structure of the singularities we restrict ourselves to the to- 
tal cross section for unpolarized W -bosons. The singular structure of the cross section for po- 
larized W-bosons is analogous.) As mentioned above the contribution involving Cq(1,2,3) 
is the only infrared divergent contribution as far as the generic loop diagrams are concerned. 
Furthermore we note that no expansion in the dimensional regulator £ has been performed so 
far. As a consequence we observe that the rational function multiplying the Cq (1,2,3) inte- 
gral is up to an additional factor the squared born matrix element in d dimensions as it must 
be. This is an important cross check of the calculation. The cancellation of the divergencies 
by the real corrections is discussed in detail in sectionHJ 

Let us now switch to the UV divergent contribution which is generated by the scalar one- 
and two-point integrals. Defining the finite parts A, B of these integrals by 

(2nv) 2£ A(m t ) = -T(l+e)(^-) m t 2 +A(m t ), 
£ \m t m\, l 



(2 » g 
— ) m b 2 +A(m b ) 
m t m}, J 



(27i//) 2e B = -T(l+e)(^-) +B , (3.12) 
£ \m t mt, J 

the UV divergent contribution before renormalization reads: 

»I*UIW = ^wg£)>„l 

- 3^r(l+e)C r (^-) -F(z„,z b ,z,), (3.13) 



2 
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with 



F(z Wl Zb,z t ) = - — I zkbt{zb-2z t ) + , - — r — (2zb (1+Zw ) -z&(l -5z w -6zw 

*Zw \ ^\ Zw)Zw 



+ 6zt{\+z w 2 ))-2 Zt {l- 5z w - 6z w 2 -2 Zt (l+ z w 2 ) ) ) 
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+ 



3(1 Zw) Zw 



(x b -X bt z) 3 (x t +X bt z) 3 

8ZbZ^ 



+ t — — rr f 64z w (zi - Zb (2zt - z w ) + z t z 



+ z t z w - 2 Zw 2 ) - T^-j^p j 10 ^ 3 +^ 2 ( 7 " 22z ' + 5 **0 - 14 ^ 2 

- lz w + 24 Zw 2 + 2 Zt (5- 7z w ) ) + (1 - z*) (2z 2 - 2z w ( 1+ z w ) - z* ( 1 - 3z w )) } 

+ 32z f w |2 Zfc 3 - z fc 2 (3 + 8z f - 3z w ) + zb{ 1 + 10z, 2 - 5z w - 2z w 2 + 2 Zt (5 - 4z w )) 

- 4z f 3 - Zt 2 (l - z w ) + z t (l - 7z w + \2 Zw 2 ) + 2z w (l + z w ) } 

- ^ + \ 8zfc Zw-12zfc (1 -z w )z w -2zi (8z, z w -Zi(3- 18z w + 31z w ) 

- 2 Zw ( 1 - 2z w - 5z w 2 ) ) - z fc (4 ( 1 - 5z w )z H 2 + 4z f 2 (3 - 27z w + 1 6z w 2 ) 
+ z t (3 - 2 l Zw + 4 lz w 2 + 25z w 3 ) ) + z, ( 8z f 3 z w + 6z, 2 ( 1 - 10z w + 9z w 2 ) 

+ 4z w (2 + 3z w - 6z w 2 + 5z w 3 ) - 3z f ( 1 - 3z w + 1 9z w 2 - 9zJ ) ) } 

- t \ — r { 16z fe 4 z w + z b 3 (6- 12z w + 66z w 2 ) -z h 2 (3-2\ Zw + 32z 2 z w + 8 lz w 2 

[Xb-AbtZ) l 

- 15z w 3 + 4z ? (3 - 30zw + 1 lzw 2 )) + Zb(z 2 (6 - 2A Zw + S2 Zw 2 ) + 4z w (l + 3z w + 4z w 3 ) 

- z f (3-21z w + 73z vl , 2 + 41z w 3 ))-8(l - Zt)z t z w (2z 2 + z w - 5z w 2 - z t (l — 3z w ))| Y 

The UV singularities must be canceled by the renormalization procedure. For the renormal- 
ization we need only to consider the wave function renormalization Zip of the quark fields 
and the renormalization of the mass parameters. As mentioned earlier the top quark mass 
and the bottom mass are renormalized in the on-shell scheme. The generic counter term for 
a quark flavour / is given by 

' X . = i((Z%(f)-l)y-(Z^(f)-l)m m ) 

= i8Z%{f)W-m m )+i(82%{f)-8Z% l {f))m m , (3.14) 

with 

Z^ )0 (/) = 1 + 8Z^ (/). (3.15) 

The first term in eq. (13.141) gives simply the corresponding born diagram multiplied by 
— 8Z^ I 1 . This contribution itself is not gauge independent, it is canceled by a similar con- 
tribution from the ytt and ybb counter terms. In addition to the mass counter term we have to 
consider the counter terms which corresponds to the vertex corrections. These counter terms 
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amount to an additional factor multiplying the born amplitude: 





W,y V = W,y / x (JZ$(f)Jz™(f)-l) 7 (3.16) 



where /(/') denote the flavours of the outgoing quark (anti-quark). Expanding the factor in 
the coupling we obtain 



Z$(f)y/Z$(f) - 1 = -(8Z5?(/) + 5Z$?(/)). (3.17) 
The contribution from the renormalization is thus given by 

Ml en = ^(dZ^(t)+8Z^(b))M 0A + ^ i ct (3.18) 

where M± denotes the contribution of the term «(8Zip — 8Zo)m on in eq. (13.141) . The contri- 
bution to the squared matrix element finally reads 

M?*Mq/ + M?"*Mq A = (SZp(0 + dZy(b))\M OA \ 2 + 2Re(M l ct Mo/). (3.19) 

Using 



e 



5Z W ) = £r(l + e)C F g) {-I + A_ 4Ko(e) 




8Z?(/)-8ZS°(/) = ^r(l+e)C f ^- {|+4^+0<£|. (3.20) 



we obtain 



ATI 



.2\ e o /A^..2\ £ O //l^..2\ e 



r(1 + e)Cf /_ W) + A« _ 8 W o/ 

I £ \m t m b J £ IR V m f 2 / £ IR V m b 2 J J 
+ ^ r (l + e)C F {~+4}(f^) f t ( Zw , Zb , Zt ) + f^\ f b (z w ,z b ,z t )) (3.21) 

In eq. (13.201 ) we have introduced £jr (with d = 4 + 2£ir) to distinguish between IR and 
UV singularities. The functions ft,f b are given in the appendix. Inspecting eq. (13.211) one 
observes that the UV singularities match exactly those appearing in eq. (I3.13I) . The UV 
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singularities are thus canceled by the renormalization procedure as it should be. Note that 
the renormalization procedure introduces an additional IR divergent contribution 



+ 




nit J V m b 

The complete IR divergent contribution of the virtual corrections is thus given by: 



a. 



r(l+e) 



s J K bt e y,r,6PoL 



^O.d| 



a. 



-^r(l + e)C F 



'4TZ/J 7 



2\ e 



+ 



4k/j 2 \ e \ 1 

m b 2 



£ l ^o,d l 



' y,t,b Pol. 



The cancellation of the IR singularities is discussed in the next section. 



4. Real corrections 



(3.22) 



(3.23) 



In this section we consider the calculation of the real corrections 

W + (k w )+y(k y ) -> t(kt) + b(k b )+g(ki). (4.1) 

The calculation of the matrix elements is straightforward and does not impose any problems. 
In principle it can also be done automatically with packages like for example CompHEP [ 53 1 
or MadGraph ll54ll . We have checked that we reproduce in the soft limit the factorization 
formulae 

M (k y ,k w ,k t ,k b ,ki) ^— > s{k t ,ki,k b ) x M (k y ,k w ,k t ,k b ) (4.2) 
with the well known eikonal factor 

2{k v kj) (krk) (kj-kj) 



(4.3) 



(ki ' ks){kj ~ k$) {ki-ks) 2 (kj-ks) 

We also compared the results numerically with MadGraph and found agreement. The IR 
divergencies arise from the phase space integration over regions where the gluon is soft. To 
extract these singularities we used the subtraction method for massive quarks ll55ll56l . The 
basic idea of the subtraction method is to add and subtract a so-called dipole term which 
on the one hand matches point-wise the singularities in the real corrections and on the other 
hand is simple enough to allow an analytic integration over the unresolved phase space in d 
dimensions 1561: 



„NLO 



dR(k t , k b ) Jo vlrt - (ky, k w , kt,k b ) + do bom (k y , k w , k t ,k b ) ® I 



Je=0 



+ dR(k t ,h,ki) do iea (ky,k w ,k t ,k b ,ki)\ £=Q 



- £ Ja born (fc Y ,fc w ,^,^)®JV dipo i e 

dipoles 



e=0. 



(4.4) 
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In eq. (14.41) the symbol ® involves in addition to the identification of the kinematics also spin 
correlations. In general also colour correlations appear. This is not the case here because the 
leading-order matrix element is proportional to the unit matrix in colour space. We note that 
the dipoles rfVdipoies are universal and can be obtained from the study of soft and collinear 
limits ll56ll . The dependence on the specific process is encoded in do bom (ky,k w ,k t ,k b ). The 
integral of the dipoles over the 'dipole phase space' which appears in the factorization of the 
phase space is denoted by 

1= f £ ^dipoles- (4.5) 
1 dipoles 

A more detailed description of the subtraction method is given in ref. [56|, here we just 
reproduce the relevant formulae for the specific case studied in this paper. In the notation of 
ref. 1 56 1 the dipoles dVdipoies are given by 

V do hom (ky,k w ,k t ,k b )®dV dipoles = 1 (Vit >b )\Mo(k lt ,k~ b )\ 2 
dixies Ah-k t ) 

+ o/, 1 , s (Vib,t)\M (kt,hb)\ 2 (4.6) 



with 



£ Vlt.fr 

l-z t {l-yit,b) nt,b 



2(h-k b ) 

2 



m t 

l+z t + - 



(k t ■ k\ 



(4.7) 



The momenta kn^kk play the role of the emitter and the spectator. For the detailed definition 
we refer to ref. Il56ll . The general expressions for Z/,y/;jfc, v,-; k, Vij.k are also given in ref. Il56ll . 
For the specific reaction considered here we obtain 

= (k ■ h) 
(kt ■ k b ) + (ki -k b y 

(kfh) 

yit ' b ~ (k r k 1 )+(h-k b ) + (k r k b y ( • } 

ht,b = — ^ , (4.10) 

l-Zt-Zb 

and 



y/ [2 Zb + (1 —zt — z b )(l —yit,b)] 2 ~^z~b , AU . 

V\t,b = y rrz ; • (4.11) 

(1 -z t -z b )(\ -y\ t ,b) 

Combining the dipole terms as given in eq. (14.41) together with the real corrections given 
by do real (kj,k w ,k t ,k b ,ki) the integration can be done numerically in 4 dimensions over the 
whole phase space. The integrals of the dipoles over the dipole phase space (which we have 
to add to account for the additional term we have subtracted from the real corrections) can 
be obtained from ref. ll56l : 

(4.12) 
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with 



Vij* 



^ln(p) -ln(p) In (l - fa + w ) 2 ) - ~]n( Pj {m,Vk)f 



|ln(p yt ( A / ; , w )) 2 + ^- + 2Li 2 (-p)-2Li 2 (l-p) 



1 2 1 2 

"2 Li2 ( 1 -PjiMhW) - 2 Ll2 ^ " PkiMhW) 

1 1 - Mk/ 



+ 0(8) 



(4.13) 



and 



Pn(Vj,Vk) 



i + ln(ju e ) - 2 In [( 1 - w ) 2 - 4] + ln( 1 - w ) 



2^ 



In 



+ 5 



2 w (l-2^ +Q(£) 



1 — A** i-M 2 Q -4 



(4.14) 



1-^^ + 2^2/(1 -1^-/4) 
\\ 1+^ + 2^/(1-^-^) 



(n = j,k) with m = y/zi 



(4.15) 



and p as defined in eq. (I3.6I) . From the formulae above we can read off the singular contri- 
bution 



a. 



1 



471/J 



2\ e 



27ir(l-e)V * J Cp e\ X bt 
Comparing the above result with eq. (13.231) we observe that the real corrections indeed cancel 
the IR divergent contribution from the virtual corrections. Having canceled the IR divergen- 
cies all the remaining phase space integrals can now be done numerically in 4 dimensions. 



2 / 1 -Zt-Zb 



ln(p) + l \Mq A 



(4.16) 



5. Results 



Before presenting the results we first discuss a few consistency checks. As mentioned earlier 
we have checked the loop integrals appearing in the virtual corrections with the FF package 
of G. J. van Oldenborgh 15UII or in the case of the box integral by comparison with re- 
sults available in the literature. Using the box integrals in d + 2 dimensions we have verified 
that only the triangle integral Cq ( 1 , 2, 3) produces an IR singularity and that the form of this 
singularity agrees with the structure predicted by QCD. By this procedure we test essentially 
the coefficients of the two box integrals and the IR divergent triangle integral. Note that this 
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Figure 5.1: Total cross section for the process W + y — ► tb for polarised W-Bosons. The 
smaller figures inside the plots show the threshold region. 



check is valid in d dimensions. That means that the y-algebra is also tested in d dimensions. 
(The treatment of 75 is not an issue here because we have only four external momenta.) Fur- 
thermore we have checked that the UV singularities have exactly the form as predicted by the 
renormalization procedure. The structure of the UV singularities is determined by the one- 
and two-point integrals. The coefficients of these integrals (more precisely a linear combi- 
nation of them) are thus checked by the fact that we reproduce the predicted structure of the 
UV singularities. We have checked the real corrections by the comparison with Madgraph. 
A further important check is also the finiteness of the real corrections in combination with 
the subtraction terms discussed in the previous section. This is a non-trivial check because 
the matrix elements are tested point wise in the singular regions. 

Let us now come to the numerical results. For the numerical evaluation we have chosen the 
following parameters: 

G f = 1.16639 x KT 5 GeV~ 2 , cc= , m w = 80.42 GeV, m t = 175 GeV. (5.1) 

137.036 

For the strong coupling we have used a running a s with the renormalization scale set to the 
center of mass energy. As input value we used 0C s (/j = 200 GeV) = 0.105. Note that the 
Fermi constant Gf and the electric coupling a enter only through a prefactor and can thus be 
changed without redoing the numerical integration. For the b-quark pole mass we consider 
the range between 4.6 and 5. 1 GeV as given by the particle data group ll57l . In fig. 15. li the 
total cross section for the process W + y — > tb is shown for polarised VK-bosons. Fig. 15 . 1 h 
shows the total cross section for longitudinally polarised W-bosons, whereas fig. 15.1b is 
the corresponding plot for the transversely polarised case. Both plots show the Born cross 
section as well as the QCD corrected cross sections for two different values of the Z?-quark 
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differential cross section 
(QCD corrections included) 
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b) transverse 

differential cross section 
(QCD corrections included) 
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Figure 5.2: Differential cross section for the process W + y — > tb with polarised W-Bosons as 
a function of cos(0 f y) 

mass. The QCD corrections are of the order of 12% in the longitudinally polarised case and 
of the order of 24% in the transversely polarised case for a center of mass energy of 1000 
GeV. In the smaller figures inside the two plots the threshold region is shown. Here one 
observes that the QCD corrections become negative for a center of mass energy between 
195 and 225 GeV. Close to the threshold the corrections are again positive. As one might 
expect from phase space arguments the cross section for a £>-quark mass of 4.6 GeV is larger 
than the cross section for a Z?-quark mass of 5.1 GeV. We note that the difference between 
the two different mass values is quite sizeable in the energy range 300-600 GeV given the 
smallness of £>-quark mass compared to the center of mass energy. The relative size in this 
region is roughly given by ln(m^ 1 /m^ 2 )/ln(m^ 1 /5). For a center of mass energy arround 
500 GeV this corresponds to an effect of around 2.2%. For larger center of mass energies 
the curves approach rapidly. Furthermore we note that the cross section for transversely 
polarised W -bosons is suppressed in comparison with the one for longitudinal W -bosons. 
This is a well known feature of the W -boson couplings to very massive quarks. In ref. II5H1 
it has been demonstrated that in the space-like axial gauge with a specific parmetrization of 
the Higgs-doublet the contribution of longitudinally polarized gauge bosons comes primarily 
from the 'scalar gauge fields'. In particular in this specific gauge the equivalence theorem 
11591 IrjUl IfjTl Ifj2l known from the i?^-gauge becomes an identity in the sense of an expansion 
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Let us add at this point also a remark about the renormalization scale uncertainty. As told 
already the QCD corrections are of leading-order in the strong coupling constant. That 
means no compensation of the residual scale dependence takes place. If we consider for 
example for a fj value of 500 GeV the range between 250 and 1000 GeV we obtain a variation 
of the QCD corrections by about 7-8 %. Keeping in mind that the QCD corrections are only 
of the order of 10-20 % this yield an uncertainty of the cross section of 1.5 %. We can thus 
conclude that the scale dependence is small as far as the total cross section is concerned. 

In fig. l5.2l the differential cross section is shown as a function of cos f y, with f y defining the 
angle between the initial state photon and the final state top quark. Figure a) shows again 
the longitudinally polarised case and figure b) the transversely polarised one. In the curves 
shown the QCD corrections are included. Both cross sections increase strongly for the case 
that the angle B t y becomes close to 1 80 degrees. The origin of this behaviour is an initial state 
collinear singularity which appears for massless ^-quarks, when the initial state photon and 
the final state £>-quark become collinear. In the case of a none-vanishing Z?-quark mass this 
singularity is regulated by the finite ^-quarks mass and becomes manifest as a ln(mj, 2 /s). 
When the photon and the top quark become collinear the increase of the cross section is 
smaller since the top quark mass is not so small compared to the center of mass energies 
considered in fig. 15.21 

In principle these logarithmic terms could be large and one might worry that the convergence 
of the perturbative expansion is spoiled. As far as QED is concerned this is here not the case 

2 

as long as one considers only moderate values for the center of mass energy. Then aln(^-) 

is still a small quantity (i.e. Tyj^ n (j^i) = —0.077) and perturbation theory remains ap- 
plicable. Although these logarithms are not a problem at leading-order it is worthwhile to 
study their resummation. This is interesting in itself for two reasons. On the one hand the 
framework to do so is the so-called structure function approach for massless quarks (only 
the Z?-quark is considered massless) — which, in principle, one could have used from the 
begining. In the structure function approach the logarithms are absorbed into the structure 
functions and resummed via an Altarelli-Parisi like evolution. As mentioned earlier the QED 
evolution itself is not important for moderate values of the center of mass energy. Therefore 
one might argue that the structure function approach for massless ^-quarks should give a 

2 

good description because terms of order ^ which are dropped in this approach — are 

small. In general this is not true because close to the threshold region the ^-quark mass ef- 
fects can be important. It is therefore instructive to compare the fixed-order calculation with 
the structure function approach. The second reason why the structure function approach is 
of interest are the logarithms appearing in the QCD corrections. In principle they could be 
larger and the need for the resummation becomes more important 1 . The theoretical frame- 
work would be once again the structure function approach but now with a mixed evolution. 

'From a practical oriented viewpoint those logarithms if present can not cause a serious problem because 
otherwise the QCD correction would be much larger. Nevertheless one should address this issue in the future 
to get a more reliable prediction. 
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In the structure function approach the total cross section for top quark production reads: 



a = J dxTy/y(fiF,x) x 6(W + (k w )y(xk y ) — > t(k t )b(k b )) 

+ J dxr b/y (/j F ,x) x a(W + (k w )b(xky) -> ?(& f )) (5.2) 

The cross sections appearing in the above equation are the subtracted cross sections for mass- 
less ^-quarks. They depend on the factorization scheme used to factorize the singular part. 
We used the MS scheme. Note that although not written down explicitly the cross sections 
depend in general also on the factorization scale /up at which the subtraction is done. The 
functions Ty/y(jiF,x), T b /y{np,x) are the 'parton distribution functions' describing the proba- 
bility to find a photon or a £-quark inside a photon. The above procedure is in fact almost the 
same as the corresponding procedure in QCD describing hadron-hadron reactions. However 
there is one important difference: in QCD the structure functions are not calculable in pertur- 
bation theory, while in QED it is possible to calculate the structure functions perturbatively. 
To the order needed here they are given by 

r y / y (ji F ,x) = 8(l-*) + 0(cc), 

r b/y (MF,x) = ^Q 2 (x 2 + (l-x) 2 )\n(^j+0(a 2 ). (5.3) 

Note that the Ty^y(jip,x) structure function is only needed to order a because the subtracted 
hard scattering cross section starts already with a. The above results for the structure func- 
tions can be easily obtained from a matching calculation. In fact one could also argue that 
one starts with the initial conditions Ty/y(m b ,x) = 8(1 —x), T b /y{m b ,x) = and generates 
the above distributions dynamically through evolution. This gives the same result. 

In the following discussion we restrict ourselves to the unpolarized cross section, the polar- 
ized case can be discussed in the same way. Using 

G(W+(k w )y(ky)^t(k t )b(k b )) = 

^l=aG f \V tb \ 2 N — -1 ((l-z r )(-2zw(ll+6z w + lUw 2 ) 



\[2 (1 — <, W ) 

(9 - 58z w + 33z w 2 - \2zJ)z t + 4(7 - 3z w 3 + 6z w 2 - 3z w )z 2 ) 
A( Zt + 2 Zw ) (1 - 2zt + Az w z t + z w 2 - Az 2 ) ln(z,) 

(l-^)V 

2 



+ (z t + 2 Zw ) (1 - 2z t - 2 Zw z t + z w 2 + 2z t 2 ) In ( v ^ ^' - ) ) ( 5.4 ) 



)) 



for the MS subtracted parton cross section we obtain the following result for the leading- 
order cross section in the structure function approach: 

a = ^-LaG f \V tb \ 2 N -^- z ((1 -z t )(-2z w (\\ + 6z w + \\z w 2 ) 
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Figure 5.3: Total born cross section for the process W + y — > tb for unpolarized W-Bosons. 
The upper figure compares the structure function approach with the fixed order calculation 
for massive ^-quarks. The lower figure shows the deviation of the structure function ap- 
proach from the massive calculation in percent. 

-(9 - 5Sz w + 33z w 2 - I2z w 3 )z t + 4(7 - 3z w + 6z w 2 - 3z w 3 )z 2 ) 
-A( Zt + 2 Zw ) ( 1 - 2 Zt + A Zw z t + z w 2 - A Zt 2 ) ln(z, ) 

+(z t + 2z w ) (\-2zt- 2 Zw z t + z w 2 + 2 Zt 2 ) In ( ^~^ ) ) (5-5) 

Note that we have used the structure functions given in eq. (15.31) — and not the evolved ones 
— which is strictly speaking only valid for /up mini,. As mentioned earlier limiting ourselves 
to center of mass energies up to 1 TeV the evolution does not change the structure functions 
very much. We have checked that the difference which one obtains using on the one hand the 
evolved structure functions and on the other hand the fixed order results (eq. &5.3\i ) is indeed 
of the order of a few per mil and thus negligible. Keeping only terms involving ln(m^) 
in the exact result eq. ( 12.201) . and dropping all terms which vanish in the limit — > 0, we 
reproduce the above result. The comparison between the two approaches is shown in fig. 15.31 
It is clearly visible that for large center of mass energies the structure function approach 
agrees with the fixed order calculation. This is due to the fact that the corrections of the 

type mi, /s can be neglected at high energies and that the logarithms of the form ln(^) are 
still of moderate size so that the resummation of these terms would not change the result. In 
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the threshold region one observes a significant difference between the two methods. In this 
region the finite b-quark mass is important because it affects the location of the threshold. 

So far we have only studied the reaction W + y — > tb which is not directly observable. While a 
high energy photon 'beam' can be realized through the interaction of low energy photons of 
a laser with high energy electrons/positrons (see for example ref. J63j) a V7-boson beam is 
not available. On the other hand one can argue that the dominant contribution to the reaction 

£ + y^tbv e (£ = e,fi) (5.6) 

proceeds via the production of an almost on-shell W-boson which then interacts with the 
photon to produce the top quark and the b-quark. This is the essence of the so called effective 
W-boson approximation J64ll65l|66l|58l|63. In this approach the intermediate W -bosons is 
considered on-shell and described through structure functions similar to the afore mentioned 
structure function approach for the ^-quark. The effective VK-approximation is similar to 
the Weizsacker- Williams [68| approximation. The total cross section for £ + y — > tb\£ in this 
approach is then given by Il64ll631l66ll58ll67l 

o(£+ y ^ f b v e ) = jdx f w + /l+ (x) a(W+y ^tb)+J dx f w + /l+ (x) a(W+y -> tb) , (5.7) 

with the structure functions fw TL given by 

oc 1 — x 



ir 



fwr/ii*) = 1 + I" ( 4) • (5-9) 

1 8 71 sin 0^, x V m vv/ 

Note that we have written down only the leading terms for the structure functions. The 'sub- 
leading' terms which are suppressed by additional powers of m w 2 /s are not universal and 
depend on the exact prescription how to define them. We note that the distribution function 
of longitudinally polarized W-bosons is very well approximated by the leading term. On 
the other hand using only the leading term for the structure function of transversly polarized 
W-bosons results in an overestimate of the cross section for energies of the order of 1 TeV. 
For the structure function fw T /e(x) we have included the sub-leading terms as given in ref. 
ll65l . An additional remark on the use of those functions is in order: while in the original 
work [65 1 a lower boundary on the allowed x values appears (x > m w /E), apparently no 
such boundary appears in refs. |Rj91 F58I. Having studied the quality of the approximation 
for center of mass energies of about 40 TeV and heavy ?,£-quarks we find that only without 
this additional constraint we obtain good agreement. For the present case we have used 
the following approach: for the contribution of the longitudinally polarized W -bosons the 
constraint x > m w /E is not used. For the transversly polarized VK-bosons we must use the 

2 

additional constraint because otherwise the distribution function (including corrections, 
see eq. (2.18) in ref. [65|) is not defined. In fig. l5.4l we show the leading order result for the 
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Figure 5.4: Total cross section for the process e + y — > tbv e . The full line shows the exact 
result while the dashed line shows the result using the effective VK-approximation. 

reaction £ + y-^ tb\£. The full line is the exact result. The exact result agrees with the result 
presented in ref. [45|. The dashed line shows the result in the effective ^-approximation 
using the above prescription. 

It is clearly visible that the accuracy of the approximation is only at the 10% level for small 
values of the center of mass energy. To obtain a reliable prediction at NLO we have com- 
bined the exact leading-order result with the QCD corrections obtained in the effective W- 
approximation. We expect that by this procedure the uncertainty due to the effective W- 
approximation is smaller than a percent and thus of the same order as the next-to-next-to- 
leading order QCD corrections. The NLO cross section obtained by this procedure is shown 
in fig. 15.51 We note that the QCD corrections which are quite sizeable at the level of the par- 
tonic reaction W + y — > tb for large values of the center of mass energy are only of the order 
of 5% for the reaction e + y — > tbv e . This is due to the convolution with the ^-distribution 
functions which gives more weight to the lower center of mass energy values. 



6. Conclusions 



In the present paper we have studied the QCD corrections for single top quark production 
in electron photon interactions. We have first calculated the QCD corrections for W + y — > 
tb. Applying the effective ^-approximation these results can be used to obtain the QCD 
corrections for e + y — > tb\ e . While the corrections are sizeable for the reaction W + y — ► tb 
they are only of the order of 5% for the reaction e + y — ► tbv e . We can thus conclude that 
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Figure 5.5: Total cross section for the process e + y — > tbv e in the effective ^-approximation. 

as far as the QCD corrections are concerned the reaction e+y — > tbv e is very well suited for 
precise measurements of the CKM matrix element V t t,. 
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A. Contributions from self-energy like counter terms 
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+ 25z w 2 + 33z w 3 ) + 16z w 4 + 12z w 2 + 4 Zw \ - — f 3 ? , 5 A 2 ^ 

J 3(1 -z w y{x t + Abtz) I 

- z b 2 (2zt - 3z w + 3)- z h (2zt 2 - 2zt (1 + 3z w ) + 5z w 2 + 2z w - 1) + 2z f 3 

- 2? 2 (Zw— 1) — Z t (z w 2 + 6z w ~ 1) +5z w 2 — Z w | +ZZb^bt 

2z b (l + 2z w ) | 128z w z fc 2 fa+z f -z w ) 
(l-z w ) 3(1 -z w ) 3 (^-X fe z) 3 
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- -ZJ-. \T7~~ — \ — Y2 1 20zb 2 z w + Zb ( 1 6z f z w - 27 z w 2 + 1 2z w - 1 ) 

J(l— ZwJ \Xb — hbtZ) < 

- 4-z t 2 z w + z t (z w 2 + Sz w — 1 ) — 2z w 3 + 2z w 2 — Az w j 

^"^^ f 2 2 2 

- tt; ttt » — z\ &Zb Zw + Zb(16ztZw + 7z w -26z w + 3)+Sz t z w 

3(1 -z w y{x b - Kbtz) L 

- z, (7z w 2 + 10z w - 1) + 10zw 3 - 1 Uh 2 + 8z w + 1 } 

- 7T7 \ i 2zb Zw ~\~ Zbi^ZtZw 2z w 2z w ) ~\~ 2Zt Zw 

3(1 -z w y{x t +M t z) I 

- Zr(zw 2 + 4 Zw -l)+2 Zw 2 }]^ +0(£ 2 ) (A.2) 
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